Fractional Drag induced in Low- Dimensional Systems by Brownian Motion of Ions in 

Liquid Flow 
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We study the frictional drag force in low-dimensional systems (2D-electron and 2D-liquid systems) 
mediated by a fluctuating electromagnetic field which originate from Brownian motion of ions in 
liquid. The analysis is focused on the [2D-system-2D-system] , [2D-system -semi-infinite liquid], and 
[2D-system-infinite liquid] configurations. We show that for 2D-electron systems the friction drag 
depends linearly on the relative velocity of the free carries of charge in the different media, but for 
2D-liquid systems the frictional drag depends nonlinear on the relative velocity. For 2D-systems 
the frictional drag force induced by liquid flow may be several orders of magnitude larger than the 
frictional drag induced by an electronic current. 

PACS: 47.61.-k, 44.40.+a, 68.35.Af 
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I. INTRODUCTION 

All media are surrounded by a fluctuating electromagnetic field because of the thermal and quantum fluctuations of 
the current density inside them. The fluctuating field is responsible for many important phenomena such as radiative 
heat transfer, the van der Waals interaction, noncontact friction and frictional drag in low-dimensional systems^. 
The relation between fluctuations and friction is determined by the fluctuation-dissipation theorem. According to 
this theorem, the fluctuating force that makes a small particle jitter will also cause friction if the particle is dragged 
through the medium. Noncontact friction and frictional drag measurements are very closely related to each other. 
In both these experiments the media are separated by a potential barrier thick enough to prevent electrons or other 
particle with a finite rest mass from tunneling across it, but allowing the interaction via the long-range electromagnetic 
field, which is always present in the gap between bodies. In noncontact friction experiments the bodies move relative 
to each other, while in frictional drag experiments the motion of charged particles (electrons or ions) is induced in 
one medium, and the reaction of free carries of charge on this motion is measured in other medium (see Fig. [1]) 

There are two mechanisms of noncontact friction and frictional drag. The electrostatic friction is due to the relative 
motion of charged bodies. This mechanism of noncontact friction was observed in^. In these experiments a charged 
probe tip oscillated close to the surface of substrate. Measurement of the damping of the probe oscillation yields 
a noncontact friction coefficient that may be related to the spectrum of the electromagnetic field fluctuation at the 
probe frequency. Using the fluctuation-dissipation theorem, the force fluctuations were interpreted^ in terms of 
near-surface fluctuating electric field interacting with static surface charge. The electrostatic friction can be related to 
the 'image' charge which is induced in the sample by the charged tip. During motion of the tip relative to the sample 
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FIG. 1: Two way for studying noncontact friction. Left: Two bodies at finite separation are sliding relative to each other. 
The noncontact friction force between bodies is mediated by a fluctuating electromagnetic field. Right: A current of charged 
particles (electrons or ions) flow in upper medium. The fluctuating electromagnetic field associated with this motion will result 
in friction force acting on the free electrons or ions in medium at the bottom. This friction force will induce an electric field 
which can be measured. 
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FIG. 2: Two bodies moving relative to each other will experience van der Waals friction due to Doppler shift of the electro- 
magnetic waves emitted by them. 
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FIG. 3: A schematic diagram of a drag experiment. A current Ji is passed through layer 1, and a voltmeter is attached to 
layer 2 to measure the induced electric field E2 due to inter layer interaction. 

the "image" charge will lag slightly behind the moving charge inducing it, and this is the origin of the electrostatic 
friction^. 

Noncontact friction exists even between neutral bodies because the fluctuating electromagnetic field originated 
from charge fluctuations in one medium will induce polarization in other medium. The interaction of the fluctuating 
electromagnetic field with the induced polarization is responsible for many important phenomena such as radiative 
heat transfer and the van der Waals interaction^. When two media are in relative motion, the induced polarization 
will lag behind the fluctuating polarization inducing it, and this gives rise to the so-called van der Waals friction. 

The origin of the van der Waals friction is closely connected with the Doppler effect. Let us consider two flat parallel 
surfaces, separated by a sufficiently wide insulator gap, which prevents particles from tunneling across it. If the charge 
carriers inside the volumes restricted by these surfaces are in relative motion (velocity v) a frictional stress will act 
between surfaces. This frictional stress is related to an asymmetry of the reflection amplitude along the direction of 
motion; see Fig. [2j If one body emits radiation, then in the rest reference frame of the second body these waves are 
Doppler shifted which will result in different reflection amplitudes. The same is true for radiation emitted by the 
second body. The exchange of "Doppler-shifted-photons" will result in momentum transfer between bodies which is 
the origin of the van der Waals friction. 

The van der Waals friction can be probed not only by measuring the friction force during relative motion of the 
two bodies, but an alternative method consists in driving an electric current in one metallic layer and studying of 
the effect of the frictional drag of the electrons in a second (parallel) metallic layer (Fig. [3]). Such experiments 
were suggested by Pogrebinskii£ and Priced, and were performed for 2D-quantum wells^^. In these experiments, 
two quantum wells are separated by a dielectric layer thick enough to prevent electrons from tunneling across it but 
allowing interlayer interaction between them. A current of density Ji = n\ev is driven through layer 1 (where n< is the 
carrier concentration per unit area in the first layer), see Fig. [3J Due to the interlayer interactions a frictional stress 
a = jv will act on the electrons in the layer 2 from layer 1 which will induce a current in layer 2. If layer 2 is an open 
circuit, an electric field E2 will develop in the layer whose influence cancels the frictional stress a between the layers. 
Experiments^ show that, at least for small separations, the frictional drag can be explained by the interaction between 



the electrons in the different layers via the fluctuating Coulomb field. However, for large inter-layer separation the 
friction drag is dominated by phonon exchange^. 

Recently, it was observed that the flow of a liquid over bundles of single-walled carbon nanotubes (SWNT) induces 
a voltage in the sample along the direction of the flo w 12 i 13 . The dependence of the voltage on the flow speed was found 
to be logarithmic over five decades of variation of the speed. There have been attempts to explain this flow-induced 
voltage in electrokinetic terms, as a result of the streaming potential that develops along the flow of an electrolyte 
through a microporous insulato r ' 15 . Earlier Krai and Shapiro 16 proposed that the liquid flow transfer momentum 
to the acoustic phonons of the nanotube, and that the resulting "phonon wind" drives an electric current in the 
nanotube. They also suggested, qualitatively, that the fluctuating Coulomb field of the ions in the liquid could drag 
directly the carriers in the nanotube. However, the first mechanism^ 6 - requires an enormous pressure^., while the 
second mechanism^ 6 - result in a very small current, of order femtoAmperes^. Ini 3 - another mechanism was proposed, 
which is related to the second idea of Ref j^ 6 -, but which requires neither localization of carriers nor drag at the same 
speed as the ions. In fact^ considered the friction between a moving point charge and the surrounding medium. This 
mechanism of friction is described by the theory of electrostatic friction 2 . For neutral systems, such as a nanotube, 
the electrostatic friction proposed ini 3 - will vanishing. 

Ir>i£ it was assumed that the liquid molecules nearest to the nanotube form a 2D-solidlike monolayer, pinned to 
the nanotube by adsorbed ions. As the liquid flows, the adsorbed solid monolayer performs stick-slip type of sliding 
motion along the nanotube. The drifting adsorbed ions will produce a voltage in the nanotube through electronic 
friction against free electrons inside the nanotube. 

Ir*i£ a model calculations of the frictional drag were presented involving a channel containing ovcrdamped Brownian 
particles. The channel was embedded in a wide chamber containing the same type of Brownian particles with drift 
velocity parallel to the channel. It was found that the flow of particles in the chamber induces a drift of the particles 
in the channel. 

In this article we study frictional drag in low-dimensional system mediated by the fluctuating electromagnetic field 
originated from Brownian motion of ions. We compare this mechanism of frictional drag with the frictional drag 
resulting from charge fluctuations in electron system. A Brief Report about this work was published ini^. 

II. VAN DER WAALS FRICTIONAL DRAG BETWEEN TWO 2D-SYSTEMS 

Let us consider two media with flat parallel surfaces at separation d -C At = ch/ksT. Assume that the free charge 
carriers in one media move with the velocity v relative to other medium. According to 20 i 21 the frictional stress between 
the two media, mediated by a fluctuating electromagnetic field, is determined by 
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where n(u>) — [exp^w/fceT — and (1 <-> 2) denotes the terms which are obtained from the preceding terms by 
permutation of indexes 1 and 2. Ri P (i = 1, 2) is the reflection amplitude for surface i for p -polarized electromagnetic 
waves. The reflection amplitude for a 2D- system is determined by 21 

Rip = y T i > ( 2 ) 
t-ip "T r 

where q p = <liriqo-i(uj,q)/ujE +1, <Ji is the longitudinal conductivity of the layer i, and e is the dielectric constant of 
the surrounding dielectric. 

For a 2D-electron system the longitudinal conductivity can be written in the form ai(uj,q) = —iojxi((jJ,q)/q 2 , where 
Xi is the finite life-time generalization of the longitudinal Lindhard response function for a 2D-electron ga o 22 i 23 

XHW ' W 1 + Wcot) X °(lo + i/r, Q)/ X °(0, Q) ' W 



where for a degenerate electron gas 



X?(w, q) = \ 2z + y/{u + z) 2 - l} , (4) 
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where r is the relaxation time, z = q/2kp, u = u)/(qvp), and kp and vf are the Fermi wave vector and Fermi velocity, 
respectively, and n s is the 2D-electron density in the layer. We define an effective electric field equal to the friction 
force per unit charge: E = cru/n s e. For v <C Vf, where vp is the Fermi velocity, the friction force depends linearly 
on velocity v. For d — 175 A at T = 3 K, and with n s = 1.5 x 10 15 m~ 2 , the electron effective mass m* = 0.067 m e , 
vp = 1.6 x 10 7 cm/s, the electron mean free path I = vpr — 1.21 x 10 5 A, and e — 10 (which corresponds to the 
condition of the experiment 9 -) we get E — 6.5 x 10~ 6 v V/m, where the velocity v is in m/s. For a current 200 nA in 
a two-dimensional layer with the width w = 20/im the drift of electrons (drift velocity v — 60m/s) creates a frictional 
electric field in the adjacent quantum well E = 4 • 10~ 4 V/m. Note that for the electron systems the frictional drag 
force decreases when the electron concentration increases. As an a example, for 2D-quantum wells with high electron 
density (n s = 10 19 m" 2 , T = 273 K, r = 4 x 10~ 14 s, e = 10, m* = m e ) at d = 175 A we get E = 1.2 x 10~ 9 u 
V/m. the thickness of the channel d c = 10 nm 

Let us now consider a fluid with the ions in a narrow channel with thickness d. For d 3> qp] 1 ^S> d c , where 
qp> = v AitNqQ 2 / £ c kBT is the Debye screening wave number (No is the concentration of ions, and e c is the dielectric 
constant of the liquid in the channel, and Q is the ion charge), the channel can be considered as two-dimensional. 
The Fourier transform of the diffusion equation for the ions (of type a) in the channel can be written in the form 

»W _ a „ 2 ( -a , N a Q 2 d c 



where u° and <j) q are the Fourier components of the surface charge density and the electric potential, respectively, and 
D a is the diffusion coefficient of the ions in the liquid in the channel. From Eq. ((5]) we get 

N a Q 2 d c 2 tp q 

a q = 1 rp g — ■ /n ■ (d) 

q k B T q 1 - iuj/D a 

The surface current density resulting from the diffusion and drift of the ions of type a, is determined by the formula 

i? = -iqD (a a q + N a Q 2 d c k B T<p q ) = 



k B T q- - iijj/D a 



where E q = —iqip q is the Fourier component of the electric field. Furthermore, there is a surface current density 
connected with the polarization of the liquid, which is determined by the formula 

e - 1 

j pq = -lUJp q = -iujdc— E q , (8) 

47T 

where p q and s c are the surface polarization and dielectric permeability of liquid in the channel, respectively. Thus 
the total current density j q = a(uj,q)E q , where the conductivity of the 2D- liquid is determined by the formula 

For the [2D-electron]-[2D-liquid] configuration with the same parameters as above for electron system (with high 
electron density) at d = 175 A, and with e c = 80, = 10 24 m~ 3 , D = 10~ 9 m 2 /s, d c = 100 A we get E = 4.6 x 10~ 8 v 
V/m, which is one order of magnitude larger than for the [2D-electron]-[2D-electron] configuration (with high electron 
density). Fig. [4] shows the dependence of the effective electric field in the channel on velocity of the liquid flow in 
ajancent channel, with the same liquid. For comparison with [2D-electron]-[2D-electron] and [2D-electron]-[2D-liquid] 
configurations, for the [2D-liquid]-[2D-liquid] configuration the effective electric field is many orders of magnitude 
larger, and depends nonlinerly on the liquid flow velocity v. 
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FIG. 4: The effective electric field in a 2D-channel with liquid as a function of the flow velocity in second 2D-channel for 
identical liquids in both channels. The temperature T — 300 K, the ion concentration in liquid No = 10 24 m -3 , the thickness of 
the channels d c — 10 nm, the diffusion coefficients of ions D = 10 _9 m 2 /s and the dielectric constant of the liquid e c = 80. The 
dielectric constant of the dielectric in the gap between channel e = 10 and the separation between the channels d = 175 A. 



III. VAN DER WAALS FRICTION AL DRAG IN A 2D- SYSTEM INDUCED BY LIQUID FLOW IN A 

SEMI-INFINITE CHAMBER 

Let us consider a 2D-electron system, isolated from a semi-infinite liquid flow by a dielectric layer with the thickness 
d. For the 2D- electron system the reflection amplitude is determined by Eqs. © - (JU). To find the reflection amplitude 
for interface between the dielectric and the liquid we will assume that the liquid fills half-space z > 0, and that the half- 
space with z < is filled by a dielectric with the dielectric constant e. Let us study the reflection of an electromagnetic 
wave from the surface of the liquid in the nonretarded limit, which formally corresponds to the limit c — * oo. In the 
region z < the potential can be written in the form 

<p q = (e-« 2 - ite^e***-***, (10) 

where q is the magnitude of the component of the wave vector parallel to surface. We will assume that the liquid 
consists of ions of two types a and b. The equation of continuity for the ions 

- ibJUi + V • J; = 0, (11) 

where i = a, b, rii = Ni — No, where and Nq are the concentration of ions in the presence and absence of the 
electric field, respectively. To linear order in the electric field 

J 4 = -NomQiVtp - AVni, (12) 

where Di is the diffusion coefficient, /ij is the mobility and Qj is the charge for ions of type i. The diffusion coefficient 
and the mobility are related with each other by the Einstein relation: Di = ksTfii. We consider the case when 
the different ion mobilities differ considerably. In this case, in the calculation of dielectric response it is possible to 
disregard the diffusion of the less mobile ions. Omitting the index i for the more mobile ions, after substitution of 
(fT2|) in ((nj we obtain 

iujn + DV 2 ( n+ ^-tp | = 0. (13) 
V k BT J 

This equation must be supplemented with Poisson's equation 

VV = 14 

where £o is the dielectric permeability of the liquid. The general solution of equations (fT3"|) and (|14p can be written 
in the form 



(15) 
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FIG. 5: The effective electric field in a 2D-channel with liquid induced by liquid flow in a semi-infinite chamber as a function 
of the flow velocity for identical liquids in the channel and in the chamber. The temperature T = 300 K, the ion concentration 
in liquid No = 10 24 m~ 3 , the thickness of the channel d c — 10 nm, the diffusion coefficients of ions D = 10 _9 m 2 /s and the 
dielectric constant of the liquid eo = 80. The dielectric constant of the dielectric in the gap between channel e = 10 and the 
separation between the channel and semi-infinite chamber d — 1 nm. 



where A = \/q 2 + qjj — ito/D and qo = \/ 'iitNoQ 2 / 'eoksT '. At the interface (z = 0) the electric potential and the 
normal component of the electric displacement field must be continuous, and the normal component of the flow density 
must vanish. From these boundary conditions we obtain 

C x + C 2 = l-R, (16) 



- e (7Ci + qC 2 ) + eq(l + R) = 0, (17) 

iuryCi + Dq 2 D qC 2 = 0. (18) 

From ([Ig ]) -([T8 ]) we get 

e- 1 

fl=-TT. ( 19 ) 



where 
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" e{Dq 2 D q-iuj\y (20) 

For v -C vf the frictional drag force acting on the electrons in the 2D- system, due to the interaction with the ions 
in the liquid, increases linearly with the fluid velocity v. In particular, for No = 10 24 m~ 3 , T = 273 K, eo = 80, 
D = 10 _9 m 2 /s, for a high electron density (n s = 10 19 m~ 2 ) in the 2D- electron system, E — 1.4 • 10 -6 t> V/m. This 
effective electric field is three orders of magnitude larger than obtained for two 2D- electron systems with high electron 
concentration, and of the same order of magnitude as friction between two 2D- electron systems with low electron 
concentration. 

Fig. \5\ shows the dependence of the effective electric field in the 2D- channel on the velocity of the liquid flow in 
the semi-infinite chamber, for identical liquid in the channel and in the chamber. We have used the same parameters 
as above for the liquid, with the separation between the channel and chamber d = lnm. The effective electric field in 
the channel initially increases with the flow velocity, reaches a maximum, and then decreases, in agreement with the 
model calculation ir>i£. The position of the maximum decreases when the density of ions decreases. The frictional 
drag force induced by the liquid flow in the narrow channel is 9 orders of magnitude larger than the frictional drag 
force induced in a 2D-electron system. 



IV. VAN DER WAALS FRICTIONAL DRAG IN A 2D- SYSTEM INDUCED BY LIQUID FLOW IN AN 

INFINITE CHAMBER 

As a limiting case of the situation considered above, let us consider a 2D-system immersed in a flowing liquid in an 
infinite chamber. We assume that the liquid flows along the n-axis, and that the plane of 2D-system coincides with 



the xy- plane. According to the fluctuation-dissipation theorem, for an infinite medium the correlation function for 
the Fourier components of the longitudinal current density is determined by^ 

< M* >«= 7A2 [™H + l/2]w 2 Ime(a; ) fc), (21) 

where k = (q, k z ) is the wave vector. The longitudinal current density is connected with the charge density pk = kjk/u; 
via the continuity equation, from which one obtain 

< PlPk >«= 7^)2 K<*>) + l/2]fc 2 Ime( W , fc). (22) 
Poisson's equation for the electric potential gives 

4 = kW^r (23) 

In the xy-plane the q - component of the electric potential is determined by 

^( W ,«,0) = y ^—A- ( 24 ) 

From ([22]) and p3j) we get 

< ^(w, ? , g, 0) >=4ft(n(o;) + l/2)ImE(a; 1 g,0), (25) 

where 

f°° rlk 1 

X(u,q,0) = - , .. . (26) 

./.^ 2vr k 2 e e(uj,k) 

Taking into account that E q — iq(f q we get < E^E*^ > UJ — q 2 < <pl<py > u - From the diffusion and Poisson's 
equations we get 

-Pk = k [ Pk+ , „ , (27) 



e k 2 ip k = 4np k + 47rp£. 



(28) 



From (|27p and (|28]) we get the dielectric function of the Debye plasma 



e{u),k) = 1 + 



k 2 - iuj/D' 



Substituting ^ in ® gives 



(29) 
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(30) 



According to the fluctuation-dissipation theorem, the average value of the correlation function for the Fourier com- 
ponents of the fluctuating surface charge density in the 2D-system is determined by^ 



f f* 

< T T 



fiq 2 



(n{uj) + l/2)Rea(Lu,q). 



(31) 



If the 2D- system is surrounded by liquid flow the electric field created by the fluctuations of the charge density 
in the fluid will induce surface charge density fluctuations in the 2D- system. The spectral correlation functions |27|) 
and pip are determined in the rest reference frame of the liquid, and of the 2D-system, respectively. In order to find 
the connection between the electric fields in different reference frames we use the Galileo transformation which leads 



to the Doppler frequency shift of the electrical field in the different reference frames. The electric field in the plane of 
the 2D- system, due to the fluctuations of the charge density in the liquid, will take the form 

E(x, t) = e-^+fl-TOt+iq-x^/ (32) 

where E q is the sum of the electric fields created by the fluctuations of the charge density in the fluid and the induced 
charge density in the 2D- system: 

E* q = E{ + iniqE(uj, q)T^tu+), (33) 

where lo + = uj + q x V and T q is the surface induced charge density. According to Ohm's law 

4 = a+El = *+(E( + 47TiqZ(u, <?)t> + )), (34) 

where a+ = a(oj + ,q) is the longitudinal conductivity for the 2D-system. The continuity equation for the surface 
charge density gives j q nd = uj + T™ d /q and from ([34]) we get 

r a <?t El 
T i _ _Z i i (35) 

q uj+ l-4iriq 2 <r+ /oj+Y,(uj,q) V ' 

and 

Ef 

q l-Amq 2 at/ui+^{ui,q)' 1 ' 

In order to find the electric field created by the charge density fluctuations in the 2D- system it is necessary to solve 
Poisson's equation in the rest reference frame of the liquid. In this reference frame the charge density takes the form 

r(x, t) = e -i("-i*V)t+i<i-* T il _ ( 37 ) 

Note that the charge density is composed from the fluctuating t' and induced r md charge density: T q * = + T q nd . 
In the presence of the liquid flow the electric field in the plane of the 2D- system, due to the fluctuating surface charge 
density, is determined by 

El 1 = A m q^-,q){r q nd +rf) 1 (38) 
where uj" = uj — q x V . From Ohm's law we get the following expression for the induced charge density 
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Substituting fl35J| in fl35J) we get 
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The friction force per unit area of the 2D-system is given by 



oo r ,2 
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where E q = E* q + E^ and t = t i + T 11 . Substituting §^ and (gDJ), gj) in gU) we get 
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FIG. 6: The effective eiectric field in a lD-eiectron system induced by liquid now in a infinite chamber, as a function of the 
flow velocity. For the same parameters for the liquid as in Fig. [5] The electron concentration per unit length in ID-system 
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3 ■ 10 m and the electron relaxation time r = 4 • 10 



Kea(uj + q x v)ImE((jj, q) 



{u + q x v) | 1 - Amq 2 o(uj + q x v)Y:{u, q)/ (uj + q x v) p 



+ (uj + q x v *-* lu) 



<l.e. f 



du[n(Lu) + 1/2] ( 



Re<r(w — q x v)ImE(u), q) 



(w - q x v) | 1 - Amq 2 (j(u - q x v)T,(uj, q)/ (w - q x v) | 5 



(cj - q x v <->• 



(43) 



where (u)±q x v <-> w) denote the terms which are obtained from the preceding terms by permutations of the arguments 
u) ± q x v and ui. With the same parameters as used above for the liquid, and for the high density 2D-electron system, 
we get E = 8.1 • 10~ 6 v V/m. For a ID-electron system we obtained a formula which is similar to Eq. (f43|) . Fig. [6] 
shows the result of the calculations of the effective electric field for a ID-electron system with the electron density 
per unit length n\ — 3 x 10 9 m _1 , the temperature T = 300 K, and with the same parameters for the liquid as used 
above. For the ID-electron system we obtained a slight deviation from the linear dependence of the frictional drag 
on the liquid flow velocity. The frictional drag for the ID-electron system is one order of magnitude larger than for 
the 2D-electron system. 

Fig. [7] shows the dependence of the effective electric field in the liquid in the 2D-channcl on the liquid flow velocity 
in the infinite chamber, assuming identical liquid in the channel and in the chamber. Qualitatively, we obtained the 
same results for a ID-channel. 

V. DISCUSSION AND CONCLUSION 

For a channel with open ends the frictional drag force will induce a drift motion of the ions in the liquid with 
the velocity Vd = D C QE /ksT '. The positive and negative ions will drift in the same direction. If ions have different 
mobility then the drifting ions will lead to an electric current whose direction will be determined by the current 
created by the ions with the largest mobility. For a channel with closed ends the frictional drag force will lead to a 
change in ion concentration along the channel. In the case of ions with the different mobilities, the friction force will 
be different for the ions with the opposite charges. As a result the ions of opposite charges will be characterized by 
different distribution functions which, as for electronic systems, will result in an electric field, and an induced voltage 
which can be measured. Let us write the friction force acting on the ions of different type in the form: F a = QE a 
and F]y = QEf,. From the condition that, in the static case, the flux density in the channel must vanish, we get 



k B T 



(<p - E a x) , 



(44) 




FIG. 7: The same as Fig. [5] but for infinite chamber. 



n b 
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These equation must be supplemented with Poisson's equation 

(ftp AttQ 
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(n a - rib). 



Substituting ([gl ft -HlS ]) in (gS]) we get 



|g =ql (2<p-AEx), 



where AE = E a — E^. The solution of (ITT)) with boundary condition 

dx x=±L/2 

where L is the channel length, has the form 
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The voltage between the ends of the channel is determined by 



U = <p(L/2) - <p(-L/2) = AE 
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(45) 
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For qr>L ^> 1 the voltage, which appears as a result of the frictional drag, will be approximately equal to U rs AEL/2. 
Furthermore, the frictional drag will induce a pressure difference Ap = nLeE. For example, if Nq = 10 24 m~ 3 , 
L = 100/im and E = lOOOV/m we get pressure difference Ap = 10 4 Pa, which should be easy to measure. Assume 
now that one type of ions are fixed (adsorbed) on the walls of the channel and an equal number of mobile ions of 
opposite sign are distributed in the liquid phase. In this case the motion of the polar liquid in the adjacent region 
will lead to frictional drag force acting on the mobile ions in the channel. For a channel with the closed ends this 
frictional drag will induce a voltage, which can be measured. 

In this paper we have shown that the van der Waals frictional drag force, induced in low-dimensional system by 
liquid flow, can be several orders of magnitude larger than the friction induced by an electron current. For narrow 
2D-and ID-channels with liquid the frictional drag force is several orders of magnitude larger than for 2D-and ID- 
electron systems. In the contrast to electron systems, the frictional drag force for a narrow channel with liquid depends 
nonlinearly on the flow velocity. These results contradict to the calculations in Ref M*, where it was assumed that the 
observed nonlinear dependence of a voltage on the liquid flow velocity was connected with the frictional drag acting 
on the electrons in the nanotubes. However, according to our calculations, the mobile ions in the channels of porous 



medium experience considerably greater frictional drag force due to van der Waals friction, than on the electrons in 
low-dimensional electronic structures. Thus, the liquid flow induced voltage observed ir*^ is more likely connected 
with the frictional drag experienced by mobile ions in an electrolyte in the channels between nanotubes, rather than 
the electrons in the nanotubes. These results should have a broad application for studying of the van der Waals 
friction and in the design of nanosensors. Such detectors would be of great interest in micromcchanical and biological 
application; ] 24 ' 25 , where local dynamical effects are intensively studied. 

Acknowledgment 

A.I.V acknowledges financial support from the Russian Foundation for Basic Research (Grant N 08-02-00141-a) 
and DFG. 



1 K. Joulain, J. P. Mulet, F. Marquier, R. Carminati, and J. J. Greffet, Surf. Sci. Rep. , 57, 59 (2005) 

2 A. I. Volokitin and B. N. J. Persson, Rev. Mod. Phys. 79 (3) (2007) 

3 B. C. Stipe, H. J. Mamin, T. D. Stowe, T. W. Kenny, and D. Rugar, Phys. Rev. Lett. 87, 096801 (2001). 

4 S. Kuehn, R. F. Loring, and J.A. Marohn, Phys. Rev. Lett. 96, 156103 (2006). 

5 A. I. Volokitin and B. N. J. Persson, Phys. Rev. Lett. 94, 086104 (2005). 

6 A. I. Volokitin, B. N. J. Persson, and H. Ueba, Phys. Rev. B 73, 165423 (2006). 

7 M. B. Pogrebinskii, Fiz.Tekh.Poluprov. 11, 637 (1977) [Sov.Phys. Semicond. 11, 372 (1977)]. 

8 P. J. Price, Physica B+C 117, 750 (1983). 

9 T. J. Gramila, J. P. Eisenstein, A. H. MacDonald, L.N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 66, 1216 (1991). 

10 U. Sivan, P. M. Solomon, and H. Shtrikman, Phys. Rev. Lett. 68, 1196 (1992). 

11 M. C. B0nsager, K. Flensberg, Ben Yu-Kuang. Hu, and A. H. Macdonald, Phys.Rev. B 57, 7085 (1998). 

12 S. Ghosh, A. K. Sood, and N. Kumar, Science 299, 1042 (2003). 

13 S. Ghosh, A. K. Sood, S. Ramaswamy, and N. Kumar, Phys.Rev. B 70, 205423 (2004). 

14 A. E. Cohen, Science, 300, 1235 (2003). 

15 S. Ghosh, A. K. Sood, and N. Kumar, Science 300, 1235 (2003). 

16 P. Krai and M. Shapiro, Phys. Rev. Lett. 86, 131 (2001). 

17 M. Das, S. Ramaswamy, A. K. Sood, and G. Ananthakrishna, Phys. Rev. E 73, 061409(2006). 

18 B. N. J. Persson, U. Tartaglino, E. Tosatti, and H. Ueba, Phys.Rev. B 69, 235410 (2004). 

19 A. I. Volokitin and B. N. J. Persson, Phys. Rev. B 77, 033413 (2008) 

20 A. I. Volokitin and B. N. J. Persson, J.Phys.: Condens. Matter 11, 345 (1999);Phys.Low-Dim.Struct.7/8,17 (1998) 

21 A. I. Volokitin and B. N. J. Persson, J.Phys.: Condens. Matter 13, 859 (2001) 

22 N. D. Mermin, Phys. Rev. B 1, 2362 (1970) 

23 F. Stern, Phys. Rev. Lett. 18, 546 (1967) 

24 R. B. M. Schasfoort, S. Schlautmann, J. Hendrikse, and A. van der Berg, Science 286, 942 (1999). 

25 N. J. Munro, K. Snow, J. A. Kant, and J. P. Landers, Clin. Chem. 45, 1906 (1999). 



